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IN THE UNITED STATES DESIGNATED/ELECTED OFFICE 
OF THE UNITED STATES PATENT AND TRADEMARK OFFICE 
UNDER THE PATENT COOPERATION TREATY-CHAPTER II 

PRELIMINARY AMENDMENT A 
PRIOR TO ACTION 

APPLICANT(S): Thomas ACKERMANN et al 

ATTORNEY DOCKET NO.: P01,0184 

INTERNATIONAL APPLICATION NO: PCT/DE99/03955 

INTERNATIONAL FILING DATE: 1 0 December 1 999 

INVENTION: METHOD AND ARRANGEMENT FOR PREDICTING 
MEASUREMENT DATA USING GIVEN 
MEASUREMENT DATA 

Assistant Commissioner for Patents, 
Washington D.C. 20231 

Sir: 

Applicants herewith amend the above-referenced PCT application, and 
request entry of the Amendment prior to examination on the United States 
Examination Phase. 

IN THE CLAIMS: 

On page 20: 

replace line 1 with -WHAT IS CLAIMED IS:-; 

Please replace original claims 1-9 with the following rewritten claims 1- 
9, referring to the mark-ups In Appendix A. 

1 . (Amended) A method for predicting measurement data until a final 
time-point using given measurement data, comprising the steps of: 

a) matching, using a processor, a stochastic process to said given 
measurement data; 

b) running simulation runs of said stochastic process from a given 
time-point until said final time-point; 

c) determining forecast measurement data for each simulation run; 

and 

d) predicting measurement data by stating a range of values, which 
is determined by said forecast measurement data, and providing said 
predicted measurement data as useable output. 
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2. (Amended) The method as claimed in claim 1, further comprising 
the steps of: 

determining a confidence range for said prediction of measurement 
data; and 

eliminating a% lowest and b% highest forecast measurement data. 

3. (Amended) The method as claimed in claim 2, wherein a% and b% 
are equal values. 

4. (Amended) The method as claimed in claim 1 , wherein said 
stochastic process is a non-homogeneous Poisson process. 

5. (Amended) The method as claimed in claim 1 , wherein said 
measurement data represents numbers of errors. 

6. (Amended) A method for predicting measurement data using given 
measurement data, comprising the steps of: 

a) matching, using a processor, a stochastic process to said given 
measurement data; 

b) sorting probability values generated by said stochastic process 
according to size, to provide a range around an expected value; and 

c) predicting measurement data within limits of said range, and 
providing said predicted measurement data as useable output 

7. (Amended) The method as claimed in claim 6, further comprising 
the steps of: 

sorting said probability values generated by said stochastic process 
symmetrically by size around said expected value. 

8. (Amended) An arrangement for predicting measurement data until a 
final time-point using given measurement data, comprising: 
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a processor unit, having a CPU, bus, memory, and input/output 
controller, configured in sucli a way that: 

a) a stochastic process can be matched to said given 
measurement data; 

b) simulation runs of the stochastic process can be carried out from 
a given time-point until the final time-point; 

c) forecast measurement data can be determined for each 
simulation run; and 

d) measurement data is predicted by stating a range of values, 
which is determined by said forecast measurement data, said measurement 
data being output in a useable form. 

9. (Amended) An arrangement for predicting measurement data using 
given measurement data, comprising: 

a processor unit, having a CPU, bus, memory, and input/output 
controller, configured in such a way that: 

a) a stochastic process can be matched to the given measurement 

data; 

b) a range can be ascertained by sorting probability values 
generated by said stochastic process according to size around an expected 
value; and 

c) said measurement data is predicted within the limits of the 
range, said measurement data being output in useable form. 

REMARKS 

The present Amendment revises the specification and claims to 
conform to United States patent practice, before examination of the present 
PCT application in the United States National Examination Phase. Pursuant 
to 37 CFR 1.125 (b), applicants have concurrently submitted a substitute 
specification, excluding the claims, and provided a marl<ed-up copy. All of the 
changes are editorial and applicant believes no new matter is added thereby. 
The amendment, addition, and/or cancellation of claims is not intended to be 
a surrender of any of the subject matter of those claims. 
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Early examination on the nnerits is respectfully requested. 



Mark Bergner 
Schiff Hardin & Waite 
Patent Department 
6600 Sears Tower 
233 South Wacker Drive 
Chicago, Illinois 60606-6473 
(312)258-5779 
Attorneys for Applicant 

CUSTOMER NUMBER 26574 



Submitted by, 
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Appendix A 
Mark Ups for Claim Amendments 

1 . (Amended^ A method for predicting measurement data until a final 

time-point using given measurement data, riF^whteh^l comprislna the steps 
of: 

3) matching, using a processor, a stochastic process [is matched 
]to ffeel said given measurement data; 

b) running simulation runs of FM said stochastic process[-afe 
eafFied-Gtrt] from a given time-point until M said final time-point; 

c) rthel determining forecast measurement data [ is dotorminod] for 
each simulation run; and 

d) predicting measurement data [ is prodiotod] by stating a range 
of values, which is determined by [thel said forecast measurement data , and 
providin g said predicted measurement data as useable output . 

2. (Amended) The method as claimed in claim 1 , \'m wh i c hl further 
comprising the steps of : 

determining a confidence range [ ie dot o rm i n o d ]for rthel said 
prediction of measurement dataf , where tho l ; and 

eliminating a% lowest and b% highest forecast measurement data[ 
are eliminated ]. 

3. (Amended) The method as claimed in claim 2, Fwl wherein 
[whieh^]a% and b% are equal values . 

4. (Amended) The method as claimed in [o ne o f the preceding claimo, 
tfwhtefelfthel claim 1. wherein said stochastic process is a non- 
homogeneous Poisson process. 

5. (Amended) The method as claimed in [ one of the preceding claimo, 
tfi-whi0fe1[the1 claim 1. wherein said measurement data represents numbers 
of errors. 



-5- 



Preliminary Amendment A 



6. IThe KAmended) A method for predicting measurement data using 
given measurement data, flfi-whtef^l comprisinq the steps of : 

a) matching, using a processor, a stochastic process [te^i=iatehed 
]to [thejsald given measurement data; 

b) [a-4Bfl§e4s~asGertaifiedr^ [tJqa~]probability values 
generated by [#ie]said stochastic process according to size, to provide a 
range around an expected value; and 

c) predicting measurement data[46-pfBd40te4] within [the~]limits of 
fthal said range , and providing said predicted measurement data as 
useable output . 

7. (Amended) T he method as claimed in claim 6, fm^ytnehl further 
comprising the steps of : 

rthel sorting said probability values generated by [thejsaid stochastic 
process[-afe^efted] symmetrically by size around [the]said expected value. 

8. (Amended) An arrangement for predicting measurement data until 
a final time-point using given measurement data, [whefeb^-l com prising: 

a processor u n it r4s-pfe¥ide4-afid1 , having a CPU, bus, memory, and 
input/output controller, configured in such a way that: 

a) a stochastic process can be matched to fthel said given 
measurement data; 

b) simulation runs of the stochastic process can be carried out from 
a given time-point until the final time-point; 

c) [the~]forecast measurement data can be determined for each 
simulation run; and 

d) measurement data is predicted by stating a range of values, 
which is determined by [thejsaid forecast measurement data , said 
measurement data being output in a useable form . 

9- (Amended) A n arrangement for predicting measurement data using 
given measurement data, rwhefebY^l comprising: 
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a processor unit[4s-pfevyed-af>41 , having a CPU, bus, memory, and 
input/out Put controller, configured in such a way that: 

a) a stochastic process can be matched to the given measurement 

data; 

b) a range can be ascertained by sorting probability values 
generated by [the]said stochastic process according to size around an 
expected value: and 

c) rthelsaid measurement data is predicted within the limits of the 
range , said measurement data being output in useable form . 
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SPECIFICATION 
TITLE 

METHOD AND ARRANGEMENT FOR PREDICTING MEASUREMENT DATA 
USING GIVEN MEASUREMENT DATA 

5 BACKGROUND OF THE INVENTION 

Field of the Invention 

[0001] The invention relates to a method and arrangement for predicting 

measurement data using given measurement data. 

Description of the Related Art 

[0002] A technical system often requires facilities for forecasting based on 
known (measurement) data, particularly in the context of error susceptibility or cost 
estimates. 

[0003] Forecasts generated by experts are generally subject to errors. 

Experts cannot carry out exact analyses, at least of highly complex systems. 

[0004] A stochastic point process, in particular a Poisson process, is 
described in Sidney I. Resnick: "Adventures in Stochastic Processes", Birkhauser 
Boston, 1992, ISBN 3-7643-3591-2, pp. 303-317 (Resnick). 

Summary of the Invention 

[0005] The object of the invention is to allow the automatic prediction 
2 0 (forecast) of measurement data using given measurement data. 

[0006] This object is achieved in accordance with the method and apparatus 
described below; developments of the invention are also described in the following 
text. 

[0007] In order to achieve this object, a method is provided for predicting 

2 5 measurement data using given measurement data, in which a stochastic process is 
matched to the given measurement data. Simulation runs are carried out from a 
given time-point until a final time-point. The forecast measurement data is 
determined for each simulation run. Measurement data for the final time-point is 
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predicted within a range of values, wliich is governed by the forecast measurement 
data. 

[0008] One development is to define a confidence range for the prediction of 
measurement data, where the a% lowest and b% highest forecast measurement 
data are eliminated. In particular, a% can equal b%. For example, a 95% confidence 
range can thus be defined by ignoring the 2.5% lowest and 2.5% highest forecast 
measurement data. 

[0009] One advantage is that the measurement data can be predicted 

(forecast) with an accuracy that is within a confidence range, from a given time- 
point This makes it possible to identify, e.g., the feasibility or impossibility of a task 
associated with the measurement data, at an early stage. Appropriate measures can 
therefore be initiated in order to counteract forecast impossibility. 

[0010] This is particularly important in the case of a complex system, e.g., a 

software development process, where the extent to which a schedule can be 
followed before the software is completed can be shown in a subsequent test phase. 
Even more important in this context is the ability to adopt countermeasures at an 
early stage if a delay has been clearly identified, e.g., in an integration test phase. 
This firstly affects the feasibility of the specified deadline (timescale) and secondly 
directly affects costs, since non-compliance with the agreed timescale often results 
in additional costs. 

[001 1] One refinement is for the stochastic process to be a non- 
homogeneous Poisson process. 

[0012] In particular, the measurement data may in one refinement comprise 

numbers of errors. This applies to software development, for example, where the 
level of maturity is documented in accordance with the errors measured in a test 
phase. Completion is directly dependent on this level of maturity. In other words, the 
software cannot be delivered to customers until most of the errors have been 
removed from the software. This is particularly important with regard to resources 
(required to test and correct errors) and costs (due to delayed delivery). 

[0013] In order to achieve the object of the invention, a method is also 

provided for predicting measurement data using given measurement data, in which 
a stochastic process is matched to the given measurement data. A range is 
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ascertained, by sorting the probability values generated by the stochastic process 
according to size, around an expected value. Measurement data is predicted on the 
basis of this range, and in particular the probability values within the range. 

[0014] One development is for the probability values generated by the 

stochastic process to be sorted symmetrically by size around the expected value. In 
particular, this means that the highest probability value represents the middle of the 
range, i.e., the expected value, whereas the next highest probability value is 
arranged to the right or left of the expected value. The next highest probability value 
is then arranged symmetrically on the other side of the expected value, in turn. 

[0015] This analytical (design) procedure provides a range, where the breadth 

of the range in turn indicates which probability values are significant in the prediction 
of the measurement data. 

[0016] In one particular refinement, the breadth of the range is determined by 

ignoring the probability values that lie below a given threshold. 

[0017] This produces a range (confidence range), which has a specific 
breadth as a result of the threshold. This breadth corresponds to the certainty with 
which the measurement data is predicted. 

[001 8] if one assumes that the stochastic process is a non-homogeneous 
Poisson process, then the non-homogeneous Poisson process defines a step size, 
particularly on a time axis t, which indicates when the next error will occur. One 
characteristic of the non-homogeneous Poisson process is that it has no memory, 
so that a "no-memory" search is carried out from each error that occurs at a specific 
time-point, for a time-point that indicates the next error. 

[001 9] In order to achieve the object of the invention, an arrangement is also 

provided for predicting measurement data using given measurement data that has a 
processor unit and is configured in such a way that: 

a) a stochastic process can be matched to the given measurement data; 

b) simulation runs can be carried out from a given time-point until a final 
time-point; 

c) the forecast measurement data can be determined for each simulation 
run; and 
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d) the prediction of measurement data for the final time-point can be 

predicted within a range of values, which is determined by the forecast 
measurement data. 

[0020] In order to achieve the object of the invention, an arrangement is 

5 further provided for predicting measurement data using given measurement data 
that has a processor unit and is configured in such a way that: 

a) a stochastic process can be matched to the given measurement data; 

b) a range can be ascertained by sorting probability values generated by 
the stochastic process according to size around an expected value; and 

10 c) the measurement data is predicted within the limits of the range. 

[0021] The arrangements are particularly suitable for carrying out the 

inventive method or the developments described above. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0022] Exemplary embodiments of the invention are shown and explained 
15 below with reference to the drawings, in which: 

Fig. 1 is a graph showing an accumulated number of errors over a test 

period; 

Fig. 2 is a graph showing the superimposed confidence ranges for different 

process models; 

2 0 Fig. 3 is a flowchart showing the steps in a method for predicting 

measurement data using given measurement data; 

Fig. 4 is a further flowchart showing the steps in a method for predicting 

measurement data using given measurement data; and 

Fig. 5 is a block diagram showing a processor unit. 

25 DETAILED DESCRIPTION OF THE INVENTION 

[0023] In order to be able to forecast a number of expected errors in a 
technical process, e.g., in a software development process, non-homogeneous 
Poisson processes (NHPP) are calibrated (i.e., matched to measurement data, such 
as the occurrence of errors over time) as follows: 

- 4 - Substitute Specification 



[0024] The following equation describes a counting process associated with 
the stochastic point process (non-homogeneous Polsson process): 



[0025] ^ 

[0026] and a time-point to defines the end of a test period, i.e., a time-point at 
which the given data ends. The stochastic processes 

^"(t)}t6R+ and (2) 

[0027] 

[0028] are searched with 

rnnom ^^^^^^ " "^^^ " ^(^0^ " ^^^H ^Cto) = hq) > a (4) , 

[0030] for all time-points where t>tO and given values a e (0.1 ) (confidence 
level) and nO e N. In particular, the following text exannines the increases in the 
stochastic countings process in relation to the time-point to. 

[0031] In the present case, where equation (1) represents a non- 

homogeneous Polsson process, the following equation (cf. Resnick) 



pKti) - H(to) . i) - .xp(- [i(ti) - i{to)]) . \M^-M. 



(5) 



[0032] 
[0033] applies for 

r^oo., 0 ^ to < ti < 00, ^ 6 No (6) 

[0034] 

[0035] and an intensity (mean measure, mean value function) of 
[0036] 

i: R"^ R"*',t i-» i(t>= EN(t) (7). 

[0037] Since the nature of the Polsson process dictates that the increases 

(error increases in this case) are Independent of previous increases, equation (5) for 

the time-points t>to to define a (minimum) range 
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[0038] 
[0039] 



[gu/ go] = [gu(tX go(t)] c Nq 

can be simplified to 



(8) 



go 



X p(N(t) - N(to) = £) ^ a I 

^=gn 

[0040] ^ 

[0041] Due to the unimodal nature of the Poisson count density, a range [gu, 

Qo] can be determined as follows: 



;9) 



10 



[0042] 

[0043] 
[0044] 

[0045] 
[0046] 

[0047] 



[0048] 

[0049] 
[0050] 

15 [0051] 



Step 1 : Sort the elementary probabilities 

: = p(N(t) - N(to) = i)r i e Nq 

into descending order and label the values sorted thus using 

P(0)' P(l)' • • • ( ^ • ® • {PO' Pi' . . .} = {p(o)' P(l)' • • •} 

P(0) ^ P(l) 2: . , . ) ; 



Step 2: 



Determine £ 



min : = 



^ e N 



0 



£ 



i = 0 



P(i) 



> a 



Step 3: Determine an index set 



I : = jio, . . . , i^min i ^ ^0 where 



PiQ t ' • • r P±£ 



min 



Substitute min{i} and go := max{i}. 

iel iel 



Step 4: 

The range from equation (8) is also referred to as the forecast range. 
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stochastic simulation (second approach) 

[0052] It is possible to determine the confidence range described using 
simulation, with the following steps: 

[0053] Step 1 : Start independent simulation runs based on the selected 

5 process model at time-point tO of the last error message m z N; 

[0054] Step 2: End a simulation run as soon as the required final time- 
point te is reached; 

[0055] Step 3: Repeat Step 2 until all simulation runs are finished; 

[0056] Step 4: Sort the numbers Ni(te) of the errors generated in the i- 

10 th simulation run in the time period (to, te), i=1, m, in descending order, and label 
the values sorted thus %)(te), . . . , N(^)(tg); and 

[0057] Step 5: Substitute 



%m.a/2j)W and 
§0 = = •^(ri..(l-a/2)l)('^e)' 



[0058] 

[0059] i.e., eliminate the (100 ■ (1-a)/2)% lowest and highest values. 

15 [0060] This produces the confidence range directly, 

[0061] Each individual simulation run is based on a simulation algorithm, 

which is known from (cf. Brately, et al., 1987): 

[0062] The simulated generation of intermediate arrival times for a non- 
homogeneous Poisson process is as follows: 

20 Stepi: Substituted := sup {X(t)}, where: 



(10) 



step 2; Generate a (pseudo) random variable X that is exponentially 
distributed with the parameter X , i.e., x := - iog(u) / X , where U is equally 
distributed over (0,1); 
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step 3: Generate a random variable U that is equally distributed over (0,1); 
and 

Step 4: If u < x{t^ + x) / X , then substitute t* := ts + X; otherwise substitute ts 
:= ts + X and go to Step 1 . 

5 [0063] The example graph in Fig. 1 shows an accumulated number of errors 
during a given test period. From time-point to, it shows a prediction range Kl for all 
time-points to + x. 

[0064] The intensity i is normally derived from equation (10) for X . For 

example the result is as follows: 

10 [0065] a) >.(t) = a ^ b • c • exp(- bt"") • t^"^ 

[0066] (;\.(t) is strictly monotonously descending for c<i, and unimodal 
for c > 1 with a definitive maximum at a point 

[0067] t_ = c^). 

[0068] b) Otherwise, X is derived in accordance with the above 

15 comments as follows: 

[0069] 

- ^ Hts\ (c < 1) V (ts > tniax) 

[0070] The graph in Fig. 2 shows the superimposed confidence ranges. In 
particular, this illustrates that possible forecasts become more scattered the further 
2 0 they extend into the future. In particular, confidence ranges calculated using 

different process models can be demonstrated in the same way as shown in Fig, 2. 

[0071] Fig. 3 shows a flowchart for the steps of a method for predicting 
measurement data using given measurement data. In Step 301, a stochastic 
process, in particular a non-homogenous Poisson process (to represent a stochastic 
2 5 count process), is matched to given measurement data. In Step 302, simulation runs 
are run from time-point to to a final time-point t© that is to be forecast. In Step 303, 
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for each simulation run, forecast measurement data is determined and a prediction 
of measurement data is restricted to a range which is covered by the measurement 
data determined by the simulation runs (see Step 304). In Step 305, a confidence 
range is determined in which a given proportion of the lowest and highest forecast 
5 measurement data is ignored in each case (this corresponds to the aforementioned 
range). The method terminates In Step 306. 

[0072] Fig. 4 shows a further flowchart for the steps of a method for 

predicting measurement data using given measurement data, in Step 401, a 
stochastic process, in particular a non-homogenous Poisson process, is matched to 
10 the given measurement data. Probability values are determined using the stochastic 
process, and these are sorted according to size around an expected value (see Step 
S 402). This sort operation results in the definition of a range, namely a confidence 

ffj range in this case. The breadth of the confidence range is determined by comparing 

El the accumulated probabilities with a given threshold. As described above, the 

hi 15 confidence range gives a distribution or uncertainty, respectively, of a time-point to in 
y the future, which allows the measurement data to be estimated in the future (see 

D Step 403). The method terminates in Step 404. 

H [0073] Fig. 5 shows a processor unit PIRZE that may be used to implement 

D the inventive method. The processor unit PRZE comprises a processor CPU, a 

^ 2 0 memory unit MEM, and an input/output interface iOS, which is used in different ways 
via an interface IFC: a graphics interface allows output to be viewed on a monitor 
MON and/or output to a printer PRT. Inputs are entered via a mouse MAS or a 
keyboard TAST. The processor unit PRZE also includes a data bus BUS, which 
provides the connection between a memory unit MEM, the processor CPU and the 

2 5 input/output interface IOS. It is also possible to connect additional components to 

the data bus BUS, e.g. additional memory, data storage (hard disk) or scanner. 

[0074] The C programming language is used in the following examples, which 
show an algorithm to define confidence ranges for forecasts and an algorithm for 
simulated definition of confidence ranges for forecasts. 

3 0 [0075] Program 1 : 

[0076] /* Definition of confidence ranges for forecasts */ 

[0077] /* based on the generalized Goel-Okomoto model */ 
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#include <stdlib.h> 
#inciude <math.h> 
#include <stdio.h> 

#deftne true 1 
#define false -1 

double mv_genGO(double,double,doub!e,double); double poisson(double,long); void kLnhppO; 
int main(argc,argv) 
int argc; 
char *argvn; 

{ 

double a.b,c,bt,stkn; 
long low,upp,lauf: 

if (argc<7) { 

pnntf{"\n\nZuwenig Argumente! \n\n"); 

printffAufruf; %s <Par1> <Par2> <Par3> <Startzeit> <End2eit>^ 
"<KNiveau>\n\n", argvfO]); return 1; 
} 

a = atof(argv[1]); 
b = atof(argv(2j); 
c = atof(argv[3J); 
bt= atof(argv[4]); 
st= atof(argvI5]); 
kn= atof(argv[6]); 

for (lauf=1;lauf< ;lauf++) { 

kLnhpp(mv_genGO,a,b,c,btbt+lauf(st-bt)/10.,kn,&low,&upp); 
printf{"Zeitpunkt: %8.2f Fehierintervall: [%d,%d]\n\ 
bt+lauf (st-bt)/10., low, upp); 
} 

return 0; 

} 

double mv_genGO{x,a,b,c) 
double x,a,b.c; 

{ retum( a*(1.0-exp(^*pow(x,c))) ); } 

double poisson{tannbda,wert) 
double lambda; 
long weft; 
{ 

long i; 

double itval,hv; 

if{lambda<600){ 
itval = exp{-iambda); 

for (i=wert;i>=1;i-) { itvai lambda/(doubte)i; } 



[0078] 
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[0079] 



else{ 

hv = exp{-lambda/(double)weft); 
itval = 1.0; 

for (i=wert;i>=1;i-) { itval *= iambda/(double)i*hv; } 
} 

return ( rtval ); 
} 

void kLnhpp(mv_nhpp, par1_nhpp, par2„nhpp, par3_nhpp, 

startjime. stop^time. k^niveau, lower, upper) double mv_nhpp(doubie,doub[e.double,doubie); 
double pari nhpp, par2 nhpp, par3_nhpp, start time, stop time, k_niveau: long •lower, *upper; 

{ 

long lauf; 

int lborder,modJow,mod_upp; 
double sum,tmp_mv, vaM, val_u; 

tmp_mv = mv_nhpp(stop_time,par1_nhpp,par2_nhpp,par3_nhpp) - 
mv_nhpp(startjime,par1_nhpp,par2lnhpp,par3_nhpp); lauf = (long)tmp_mv; 
*lower = lauf-1 ; 

*upper = lauf+1; modJow= false; mod_upp= false; sum poisson(tmp_mv,lauf); valj = 

poisson(tmp_mv.*lower); vaLu = poisson(tmp„mv,*upper); 
white (sum<k_niveau) { 
if {valj<val_u) { 
sum += val_u; 
(*upper)++; 
Iborder = false; 
mod^upp = true; 

vaLu = paisson{tmp mv,*upper); 

} 

else { 

sum += valJ; 
(*lower)-; 
Iborder - true; 
mod Jow = true; 

vaij = poisson(tmp mv,*lower): 

} 

> 

if (Iborder == true) { (*lower)++; } 
else { Cupper)-; } 

if (mod Jow == false) { {*lower)++; } 
if (mod_upp =- false) { (*upper)-; } 

return; 

} 

Program 2: 



[0080] 



/* Simulated definition of confidence ranges for forecasts */ 



[0081] 



/* based on the generalized Goel-Okomoto model 7 



- 11 - 



Substitute SPECiFiCATiON 



#include <stc(lib.h> 
#inciude <math.h> 
#include <time,h> 
#include <stdio.h> 
#includ6 <vaiues.h> 

#define taie 1 
#define faise -1 

double drand48(votd); 
void srand48(iong); 

double sim_exp(dQubie); double lambda^enGO{doubie,double,double,doubie); void stm_nhpp(); 

int matn(argc,argv) 

int argc; 

char *argv[]; 

{ 

timejt; double a,b.c,btst,pnt[1000000],check time(12]: longlauf,no _pnt,seed_run; intclauf; 
FILE *datei; 
if (argc<6) { 

pnntf{"\n\n2uwenig Argumente! \n\n"); 

printf('*Aufruf: %s <Paf1> <Par2> <Par3> <Start2eit> <End2eit>\n\n", 
argv[0]); return 1; 
) 

date! = fopen("sim.seed","r"); 
if (dat8i==NULL) { 
seed^run = 1; 
} 

else { 

fsca^f(datel,"%6d^&seed_^un); fclose(datei); seed run++; 
} 

datei = fopenf sim.seed^''w+''); 
fprintf(datei, "%6d\n", seed.run ); 
fclase(datei); 

time (&t) ; r tnttialisierung des */ 

t += seed_run*100 ; /* Zufallszahlengenerators V 
srand48 ((unsigned long) t) ; /* mit Hilfe der Systernzeit */ 
a = atof(argv[1]); 
b = atof(argv[2I); 
c = atof(argv[31); 
bt= atof(argv[4I); 
st= atof(argv[5]); 

sim_nhpp(lannbda^enGO,a,b,c,bt.st,&pnt,&no_pnt); 
for (lauf=1;lauf<=no J3nt;!auf++) { 
pnntf("%i5.7f %iOd \n^ pnt[iauf|. lauf); 
) 



[0082] 
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[0083] 

datei « fopen("ki.tmp',*a*); 
for (Iauf=1;lauf< ;lauf++) { 
check_time[laufl = bt+lauf (st-btyiO.; 
} 

check_time[1 1) = pnt(no_pnt]+1; r grOOer ats die grente 
simuHerte Zeit V 
ctauf = 1 ; 

for (iauf=1;iauf<=no _pnt;lauf++) { 

while (pntllauf]>=check_time[claufl) { fprintf{d3tei, "%8.2f %6d check_tfme[daafl, fauf-1); clauf++; 
} 

} 

if (pnt[no ^nt] <check_tlmel101) { 

for (Iauf=clau1:lauf< ;lauf++) { fpnntf(datej, "%8.2f %6d cfieckjimellauf], no_pnt); 
} 

} 

fprintf(datei. 'W); 
fciose(datei); 

return 0; 
} 

double sim_exp(lambda) 
double lambda; 

{ return( -log(drand48())/iambda ); } 

double iambda_jgenGO{x,a,b,c) 
double x,a,b,c; 

{ retum( a*b*c*pow(x.c-1)*exp(-b*pow{x,c)) ); } 

void sim_nhpp(lambda_nhpp, par1_nhpp, par2_nhpp, par3_nhpp, 

start_time, stopjme, path, no_j)oints) double lambda_nhpp(double,double.double,double); double 
parl_nhpp, par2 nhpp, par3 nhpp, start time, stop time; double pathQ; long *no_j>otnts; 

{ 

double sim_time,x»u,x_bar,tambda_bar; 

*no_points=0; 

stm_tim6 = start_time; 

do{ 

if (par3_nhpp<=1) { lambda bar = lambda_nhpp(sim time,par1 nhpp,par2_nhpp,par3 nhpp); 
} 

else { 

x_bar = pow((par3_nhpp.1.0)/par2_nhpp/paf3_nhpp,1.0/par3_nhpp); 
if (sim Jime>=x_bar) { " 
iambda.bar « lambda_nhpp(s'tm time,par1 nhpp,par2_nhpp,par3 nhpp); 
} 

else { 

lambda_bar - lambda.nhpp(x bar,par1 nhpp.par2_nhpp,par3_nhpp); 
} 

} 

X - sim_exp{Iambda_bar); 
u = drand48(); ~ 



if (u<-lambda_nhpp(sim_time+x,pan_nhpp»par2„nhpp,par3_nhpp)/lambda_bar) { (*no_points)++; 
pathCnojoints]=sim_time+x; 

} 

stm_time+=x; 
} 

while (sim_time<-stop_tinfie); 

return; 

} 
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[0084] 
[0085] 
[0086] 



Program 3: 

/* Definition of confidence ranges from the simulation data */ 
/* (the simulation data is sorted into ascending order) */ 



O 

m 
m 
m 
m 
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#include <stdlib,h> 
#inciude <math,h> 
#include <stdio.h> 

int qsortJcmp(int*,int*); 
int qsortJcmp(x,y) 
int "x, *y; 
{ 

if (*x<*y) { return ( -1 ); } 
else if {*x==*y) { return ( 0 }; } 
else { return ( 1 ); } 

} 

int main(argc,argv) 
int argc; 
char *argv[]; 

{ 

int pnt[11][1000QO]; 
int qsI1 00000]; 
char *dname; 
int frac,i; 

long iauf,iower_bound,upper_bound; 
long l,nojnt,seed_run; 
double ctinne[11],x: 
FiLE *datei; 



if (argc<3) { 

printf("\n\nZuwenig Argumente! \n\n"); pnntf("Aufruf: %s <Dateiname> <Konfidenzniveau (in 
%%)>\n\n-. argv[0]); return 1; 

} 

dname = argv[1]; 

frac = 100-atoi{argv[2]); 

lauf^O; 

date! = fopen(dname, V); 
if (datei==NULL){ return 1;} 
else{ 

while (Ifeof(datei)) { 
lauf++; 

for (i=:1;i<=9;i++}{ 

f5canf(datei/%8if%6d". &ctinrje[i], &pnt[i][iauf]); 

fscanf(datei/%8lf %6d\n-, &ctime[10], &pntt10][lauf]): 

fciose(datei); 
} 

lower_bound = (Iong)fioor(lauf*frac/200.); upper_bound = (Iong)ceil(lauf*(200.-frac)/200.); 
if (lower_bound<1) {lower_bound = 1;} 

[0087] 
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printf("\n\n%2d%%-Slcherheitsbereich bei %d Simulationslaufen\n\n". 

100-frac,lauf); 

for (f=i;i< ;!++){ 

for (}=1;l<=lauf:l++){ 
qs[\] = pntti][l]: 
} 

qsort(&qs[1], lauf. sizeof(lnt), &qsortJcmp); 
printfC'Zeitpunkt: %8.2f Fehlerintervall: [%d,%d]\n". 
ctime[i], qs[lower_bound], qs[upper_boundI); 
} 

return 0 ; 
} 



[0088] The above-described method and apparatus are illustrative of the 
principles of the present invention. Numerous modifications and adaptations will be 
readily apparent to those skilled in this art without departing from the spirit and 
scope of the present invention. 
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ABSTRACT 

[0089] A method and an arrangement are provided for predicting 
measurement data using given measurement data, in which a stochastic process is 
matched to the given measurement data. Simulation runs are carried out from a 
given time-point until a final time-point. The forecast measurement data is 
determined for each simulation run. Measurement data for the final time-point is 
predicted within a range of values, which is determined by the forecast 
measurement data. 
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OV/ 868239 

JC18Rec*dTOPT0 1 Sjm 2005 

rOescriptionl SPECIFICATION 

[ M e t ho d a nd a r rangem e nt for pred i ct i ng measurem e nt data usin g given 

measurement data ] 

TITLE 

5 METHOD AND ARRANGEMENT FOR PREDICTING MEASUREMENT DATA 

USING GIVEN MEASUREMENT DATA 

BACKGROUND OF THE INVENTION 

Field of the Invention 

rOOOn The invention relates to a method and arrangement for predicting 
10 measurement data using given measurement data. 

Description of the Related Art 

r00021 A technical system often requires facilities for forecasting based on 
known (measurement) data, particularly in the context of error susceptibility or cost 
estimates, 

15 r00031 Forecasts generated by experts are generally subject to errors. 
Experts cannot carry out exact analyses, at least of highly complex systems. 

r00041 A stochastic point process, in particular a Poisson process, is 
described in rf441 Sidnev 1. Resnick: "Adventures4 fh in Stochastic Processes", 
Birkhauser Boston, 1992, ISBN 3-7643-3591 pp. 303-317 (Resnick) . 

2 0 Summary of the Invention 

rOOOSI The object of the invention is to allow the automatic prediction 
(forecast) of measurement data using given measurement data. 

r00061 This object is achieved in accordance with the [features of the 
ind e p e nd en t patent claim s. Devel opmentsl method and apparatus described 
2 5 below: developments of the invention are also described in the 
f d e p e ndent l following r claimsl text . 

r00071 In order to achieve [tl^]th|s object, a method is rproposedl provided 

for predicting measurement data using given measurement data, in which a 
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stochastic process is matched to the given measurement data. Simulation runs are 
carried out from a given time-point until a final time-point. The forecast 
measurement data is determined for each simulation ain. Measurement data for the 
final time-point is predicted within a range of values, which is governed by the 
5 forecast measurement data. 

10008] One development is to define a confidence range for the prediction of 
measurement data, where the a% lowest and b% highest forecast measurement 
data are eliminated. In particular, a% can equal b%. For example, a 95% confidence 
range can thus be defined by ignoring the 2.5% lowest and 2.5% highest forecast 
10 measurement data. 

m 100091 One advantage is that the measurement data can be predicted 

m (forecast) with an accuracy that is within a confidence range, from a given 

time-point . This makes it possible to identifv- e.Q v. e.g.. the feasibliitv or 
W Impossibi lity of a task associated with the measurement data, at an early 

r 15 stage. Ap propriate measures can therefore be initiated in order to counteract 
y forecast impossibility. 

100101 [ One advantago is that tho m o asurem e nt dat a can bo prodictod 
P (for e cast) with an accuracy that is within a confidonco range, from a giv e n time 

po i nt . T h is m a kes i t pos s ib l e to idontify ][-evgT][ t he fe as i bility or impossibility of a taok 
2 0 associat o d w it h th o m e as u rem o n t data, a t an early stag e . App r opria te m e asures can 
therefore be initiated in ord e r to count e ract forecast impossibility. ]This is particularly 
important in the case of a complex system, e.g.^ a software development process, 
where the extent to which a schedule can be followed before the software is 
completed can be shown in a subsequent test phase. Even more important in this 

2 5 context is the ability to adopt countermeasures at an early stage if a delay has been 

clearly identified, e.g.^ in an integration test phase. This firstly affects the feasibility 
of the specified deadline (timescale) and secondly directly affects costs, since non- 
compliance with the agreed timescale often results in additional costs. 

room One refinement is for the stochastic process to be a non- 

3 0 homogeneous Poisson process. 

100121 In particular, the measurement data may in one refinement comprise 
numbers of errors. This applies to software development, for example, where the 
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level of maturity is documented in accordance with the errors measured in a test 
phase. Completion is directly dependent on this level of maturity. In other words, the 
software cannot be delivered to customers until most of the errors have been 
removed from the software. This is particularly important with regard to resources 
5 (required to test and correct errors) and costs (due to delayed delivery). 

F00131 In order to achieve the object of the invention , a method is also 

f proposed l provided for predicting measurement data using given measurement 
data, in which a stochastic process is matched to the given measurement data. A 
range is ascertained, by sorting the probability values generated by the stochastic 
10 process according to size, around an expected value. Measurement data is 

predicted on the basis of this range, and in particular the probability values within the 
range. 

r0014l One development is for the probability values generated by the 
stochastic process to be sorted symmetrically by size around the expected value. In 
15 particular, this means that the highest probability value represents the middle of the 
range, i.e.^ the expected value, whereas the next highest probability value is 
arranged to the right or left of the expected value. The next highest probability value 
is then arranged symmetrically on the other side of the expected value, in turn. 

r00151 This analytical (design) procedure provides a range, where the breadth 
2 0 of the range in turn indicates which probability values are significant in the prediction 
of the measurement data. 

r00161 In one particular refinement, the breadth of the range is determined by 
ignoring the probability values that lie below a given threshold. 

r00171 This produces a range (confidence range), which has a specific 

2 5 breadth as a result of the threshold. This breadth corresponds to the certainty with 

which the measurement data is predicted. 

r00181 if one assumes that the stochastic process is a non-homogeneous 
Poisson process, then the non-homogeneous Poisson process defines a step size, 
particularly on a time axis t, which indicates when the next error will occur. One 

3 0 characteristic of the non-homogeneous Poisson process is that it has no memory, 

so that a "no-memory" search is carried out from each error that occurs at a specific 
time-point, for a time-point that indicates the next error. 
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r00191 In order to achieve the objec t of the invention , an arrangement is also 
[ p r G po s ecn provided for predicting measurennent data using given measurement 
datafv w here by! that has a processor unit [is provided land is configured in such a 
way that: 

5 a) a stochastic process can be matched to the given measurement data; 

b) simulation runs can be carried out from a given time-point until a final 
time-point; 

c) the forecast measurement data can be determined for each simulation 
run: and 

10 d) the prediction of measurement data for the final time-point can be 

predicted within a range of values, which is determined by the forecast 
measurement data. 

r00201 in order to achieve the object of the invention , an arrangement is 
further [pfepeee^lEroyided for predicting measurement data using given 
15 measurement dataf r w he r ebvl that has a processor unit [is provlde4-]and Js 
configured in such a way that: 

a) a stochastic process can be matched to the given measurement data; 

b) a range can be ascertained by sorting probability values generated by 
the stochastic process according to size around an expected value: and 

2 0 c) the measurement data is predicted within the limits of the range. 

r002n The arrangements are particularly suitable for carrying out the 
inventive method or the developments described above. 

BRIEF DESCRIPTION OF THE DRAWINGS 

f00221 Exemplary embodiments of the invention are shown and explained 
25 below with reference to the drawings, in which: 

Fig. 1 is a graph showing an accumulated number of errors over a test 

period; 

Fig. 2 is a graph showing the superimposed confidence ranges for different 

process models; 
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Fig. 3 is a flowchart showing the steps in a method for predicting 

measurement data using given measurement data; 

Fig, 4 is a further flowchart showing the steps in a method for predicting 

measurement data using given measurement data: and 

5 Fig. 5 rsbowsl is a block diagram showing a processor unit. 

DETAILED DESCRIPTION OF THE INVENTION 

r00231 In order to be able to forecast a number of [ errors to be ]expected 
errors in a technical process, e.g.^ in a software development process, non- 
homogeneous Poisson processes (NHPP) are calibrated^ (i-e.^ matched to 
10 measurement data, Fe^l such as the occurrence of errors over timeHl as follows: 

r00241 The following equation describes a counting process associated with 
the stochastic point process (non-homogeneous Poisson process): 

roo25i 

r00261 and a time-point to defines the end of a test period, i.e.^ a time-point at 
15 whicJi the given data ends. The stochastic processes 

{U(t)}^^jj+ and <2) 
{0{t)}t^R+ (3) 

100271 

r00281 are searched with 

p(u(t) < N(t) - N(to) ^ 0(t)j N(to) = no) S a (4) , 

r00291 ^ 

r00301 for all time-points where t>tO and given values a e (0.1 ) (confidence 
2 0 level) and nO e N. In particular, the following text examines the increases in the 
stochastic countings process in relation to the time-point to. 

rooan in the present case, where equation (1 ) represents a non- 
homogeneous Poisson process, the following equation (cf. FWI Resnick ) 



-5- 



Mark-Up forSubstitute Specification 



p(N(ti) - N(to) . = exp(- [l(ti) - i(to)]) . 



(5) 



r00321 

r00331 applies for 

0 < to < ti < 00, ^ e No (6) 

r00341 

r00351 and an intensity (mean measure, mean value function) of 
r00361 

i: R-^ R+,t H)' i(t)= EN(t) (7) . 

r00371 Since the nature of the Poisson process dictates that the increases 
(error increases in this case) are independent of previous Increases, equation (5) for 
the time-points t>to to define a (minimum) range 

[gu.go] - [gu(tXgo(t)] c Nq (S) 

r00381 

r00391 can be simplified to 

go 

£p(N(t) - N(to) = ^) ^ a (9) . 

^=gu 

mm 

r00411 Due to the unimodal nature of the Poisson count density, a range [gu, 
go] can be determined as follows: 

r00421 step 1 : Sort the elementary probabilities 

: = p(N(t) - N(to) = (), i e Nq 

£00431 

r00441 into descending order and label the values sorted thus using 

P(0)' P(l)' • • • {PO' Pi' • • •} = {p(0)' P(l)' ■ • •} 

P(0) ^ P(l) ^ • • • ) ' 
r00451 ^ ' 
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r00461 
r00471 



step 2: 
Step 3: 



Determine 

Determine an index set 



Z P(i) ^ " 
i = 0 



r00481 
r00491 

5 roosol 
roosn 



I : = jio, . . . , i^rnin} ^0 where 



r • • • r 



Pi^ 



mxn 



P(0)' 



^(^min) J 



Step 4: 



Substitute := min{i} and := Tnax{i} 



iel 



The range from equation (8) is also referred to as the forecast range. 



Stochastic sinnulation (second approach) 

F00521 It is possible to deternnine the confidence range described using 
simulation, with the following steps: 

10 r00531 Step 1: Start independent simulation runs based on the selected 

process model at time-point tO of the last error message m e N; 

r00541 Step 2: End a simulation run as soon as the required final time- 
point t© is reached; 

[00551 Step 3: Repeat Step 2 until all simulation runs are finished; 

15 r00561 step 4: Sort the numbers Ni(te) of the errors generated in the i- 
th simulation run in the time period (to, te), i=1, m, in descending order, and label 
the values sorted thus [7]%)(te)/ . . w N(^)(t3)i_and 

r00571 step 5: Substitute 



r00581 
20 r00591 



\m-a/2j)W and 
5° = = "(rm-(l-a/2)l)(^e). 



i.e.i eliminate the dOOf^ l- n-aV2)% lowest and highest values. 
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r00601 This produces the confidence range directly. 

r006n Each individual simulation run is based on a simulation algorithm, 
which is known from (cf. [MI Bratelv, et aL, 1987 ): 

r00621 f Si m ulatod l The simulated generation of intermediate arrival times for 

5 a non-homogeneous Poisson process is as follows : 

Step1: Substitute a; := sup {X(t)}, where: 

t>tc. 



(10) 



X(t):=§i 
dt 

Step 2: Generate a (pseudo) random variable X that is exponentially 
distributed with the parameter X, \.e.^ x := - log(u) / X , where U is equally 
10 distributed over (0,1 

Step 3: Generate a random variable U that is equally distributed over (0,1 )[T]i 
and 

Step 4: If u < x{t^ + x)/X, then substitute t* := ts + X; otherwise substitute ts 
:= ts + X and go to Step 1 . 

15 f00631 The example graph in Fig. 1 shows an accumulated number of errors 
during a given test period. From time-point to, it shows a prediction range Kl for all 
time-points to + x. 

r00641 The intensity i is normally derived from equation (10) for X . For 
example the result is as follows: 

20 r00651 a) ;.(t) = a . b . c . exp(- bt^) • t^~^ 

[00661 (A,(t) is strictly monotonously descending for [ , a nd un imodal for ] [ 
w i th a d e f in i tive ma xim u m at a po i nt ] c < i , and unimodal for c > i with a 
definitive maximum at a point 



100671 t_ = c^^). 

25 r00681 b) Otherwise, X is derived in accordance with the above 

comments as follows: 
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X{ts)r (c < 1) V (ts > tjiiax) 




roo7oi 



The graph in Fig. 2 shows the superimposed confidence ranges. In 



10 



2 0 



25 



particular, this illustrates that possible forecasts become more scattered the further 
they extend into the future. In particular, confidence ranges calculated using 
different process models can be demonstrated in the same way as shown in Fig. 2. 



measurement data using given measurement data. In Step 301 , a stochastic 
process, in particular a non-homogenous Poisson process (to represent a stochastic 
count process), is matched to given measurement data. In Step 302, simulation runs 
are run from time-point to to a final time-point te that is to be forecast. In Step 303, 
for each simulation run, forecast measurement data is determined and a prediction 
of measurement data is restricted to a range which is covered by the measurement 
data determined by the simulation runs (see Step 304). In Step 305, a confidence 
range is determined in which a given proportion of the lowest and highest forecast 
measurement data is ignored in each case (this corresponds to the aforementioned 
range). The method terminates in Step 306. 

r00721 Fig, 4 shows a further flowchart for the steps of a method for 

predicting measurement data using given measurement data. In Step 401, a 
stochastic process, in particular a non-homogenous Poisson process, is matched to 
the given measurement data. Probability values are determined using the stochastic 
process, and these are sorted according to size around an expected value (see Step 
402). This sort operation results in the definition of a range, namely a confidence 
range in this case. The breadth of the confidence range is determined by comparing 
the accumulated probabilities with a given threshold. As described above, the 
confidence range gives a distribution or uncertainty, respectively, of a time-point to in 
the future, which allows the measurement data to be estimated in the future (see 
Step 403). The method terminates in Step 404. 



roo7n 



Fig. 3 shows a flowchart for the steps of a method for predicting 
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r00731 Fig. 5 shows a processor unit PRZE that may be used to implement 
the inventive method . The processor unit PRZE comprises a processor CPU, a 
memory unit MEM, and an input/output interface lOS, which is used in different ways 
via an interface IPC: a graphics interface allows output to be viewed on a monitor 
5 MON and/or output to a printer PRT. Inputs are entered via a mouse MAS or a 
keyboard TAST. The processor unit PRZE also includes a data bus BUS, which 
provides the connection between a memory unit MEM, the processor CPU and the 
input/output interface lOS. It is also possible to connect additional components to 
the data bus BUS, e.g. additional memory, data storage (hard disk) or scanner. 

10 [00741 The C programming language is used in the following examples, which 
show an algorithm to define confidence ranges for forecasts and an algorithm for 
simulated definition of confidence ranges for forecasts. 

r00751 Program 1: 

r00761 /* Definition of confidence ranges for forecasts */ 
15 r00771 /* based on the generalized Goel-Okomoto model */ 
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#inciude <stdlib.h> 
#inc!ude <math.h> 
include <stdio.h> 

#define true 1 
#define false -1 

double mv_genGO(doubie,doubie.double,double); double poisson(doub[e,tong); void ki_nhpp(}; 

int main(argc,argv) 

int argc; 

char *argv[]; 

{ 

doublea.b,c,btst,kn; 
long low.uppjauf; 

if (argc<7) { 

printf("\n\nZuwenig Argumente! \n\n"); 

printf("Aufruf: %s <Par1> <Par2> <Par3> <Startzeit> <Endzeit>", 
*'<KNiveau>\n\n\ argvfOJ); return 1; 
} 

a = atof{argvI1]); 
b = atof(argv[2]); 
c = atof(ai^v[3]); 
bt= atof(argv[4]); 
st= atof(argv[51); 
kn= atof(argv[6]); 

for (iauf=1;lauf< ;Iauf++) { 

ki_nhpp(mv_genGO,a,b,c,btbt+laur(st-bt)/10.,kn,&low,&upp): 
printf("Zeitpunkt: %8.2f Fehterintervall: {%d.%d]\n^ 
bt+iauf (st-bt)/10., low. upp); 
} 

return 0; 

} 

double mv_jgenGO{x,a,b,c) 
double x.a.b.c; 

{ retum( a*(1.0-exp(-b*pow(x,c))) ); } 

double poisson{iambda,w€rt) 
double tambda; 
long wert; 
{ 

long i; 

double rtval.hv; 

if (lambda<600) { 
itvai = exp(-Iambda); 

for (i=wert;i>=1;H { itvai *= iambda/(double)i; } 
} 

r00781 
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r00791 



e!se{ 

hv = exp(-iambcla/(double)wert); 
itval = 1.0; 

for (i=wert;i>=1;i-) { itval *= iambday(doub!e)i*hv; } 
} 

return { itval ); 
} 

void ki_nhpp(mv_nhpp, par1_nhpp, par2_nhpp, par3_nhpp, 

start_time. stop_time, k^niveau, iower, upper) double mv_nhpp(double,double,double,doubie); 
double pan_nhpp, par2 nhpp, par3 nhpp. start time, stop.time, k_niveau; long *lower, *upper; 

{ 

long lauf; 

int lborder,modJow,mod_upp; 
double sum,tmp_mv, vaij, vaLu; 

tmp_mv = mv_nhpp(stop_time,par1_nhpp,par2_nhpp,par3_nhpp) - 
mv_nhpp(start_time,par1_nhpp,par2_nhpp,par3_nhpp); lauf = (long)tmp_mv; 
*lower = lauf-1 ; 

*upper = lauf+1; modJow= false; mod_upp= false; sum = potsson(tmp_mvJauf); valJ = 

poisson(tmp_mv,*lower); vaLu = poisson(tmp_mv,*upper); 
while (sum<k_niveau) { 
if (valj<vai_u) { 
sum vai^u; 
(*upper)++; 
(border - false; 
mod^upp - true; 

val_u = poisson(tmp_mv,*upper); 
} 

else { 

sum valj; 
(*lower)-; 
Iborder - true; 
mod Jow = true; 

valj = poisson(tmp_mv,*{ower); 

} 

} 

if (Iborder == true) { (*lower)+-»-; } 
else { Cupper)-; } 

If (mod Jow = false) { (*lower)++; } 
if (mod^upp ~ false) { (*upper)-: } 

retum; 

} 

Program 2: 

roosOl /* Simulated definition of confidence ranges for forecasts */ 
roosn /* based on the generalized Goel-Okomoto model */ 
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^^nclude <stdlib.h> 
#indude <math.h> 
#include <time.h> 
#include <stdio.h> 
#tnclude <vaiues.h> 

#define true 1 
#define false -1 

double drand48(void); 
void srand48(iong); 

double sim_exp(doubie); double lambda_^enGO(double,double,double»double); void sim_nhpp{); 

int main(argc,argv) 

int argc; 

char •argvQ; 

{ 

time J t; double a,b,c,btst,pnt(1000000],check_time[12]; long lauf.no _pntseed_run; int clauft 
FILE *datei; 
if (argc<6) { 

printf("\n\n2uv/enlg Argumentel \n\n"); 

printf("Aufruf: %s <Par1> <Par2> <Par3> <Startzeit> <End2eit>\n\n", 
argv[0]); return 1; 



datei = fopen{"sim.seed",**r}; 
if {datei==NULL) { 
seed^run = 1 ; 
} 

e[5e{ 

fscanfCdatei/yoSd-.&seed^run); fclose(datei); seed_run++; 
} 

datei = fopen("S!m.seed",''w+''); 
fprintf{datei. "%6d\n", seed_run ); 
fclose(datei); 

time (&t) ; r Initialisierung des 7 

t += seed_run*100 ; /* Zufaliszahlengenerators */ 

srand48 ({unsigned long) t) ; /* mit Hilfe der Systemzeit V 

a = atof(argv[1]); 

b = atof(argv[2I); 

c = atof(argv[31); 

bt= atof(argv[4]); 

st= atof(argv[51); 

sinn_nhpp(!ambda_genGO,a,b,c,btst&pnt,&no_pnt); 
for (laijf=;1;lauf<=no _pnt;iauf++) { 
printf("%15.7f %lOd \n", pnt[iaufl, lauf); 
} 



r00821 
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r00831 

datei = fopen("ki.tmp",''a"); 
for (lauf=1 ;lauf< ;lauf++) { 
check_time[lauf] = bt+]aur(st-bt)/10.; 

} 

check_time[1 1] = pnt[no_pnt]+1 : /* groOer als die grOEte 
simuiierte Zert */ 
clauf = 1 ; 

for (lauf=1;iauf<=no j3nt;iauf+-»-) { 

while (pnt[laufl>=:check_time[claufj) { fprintf(datei. "%8.2f %6d checMmeldauf]. lauf-1); dauf++: 
} 

if (pnt[no^nt] <check_timeI10]} { 

for {lauf=clauf:lauf< ;lauf++) { fpnnt^datei. "%8.2f %6d checkjime[iaufj. no _pnt); 
} 

fprlntf(datei. "\n"); 
fclose(datei); 

return 0; 
} 

double sim_exp(lannbda) 
double lambda; 

{ retum( -log(drand48())/lambda ); } 

double lambda_genGO(x,a,b,c) 
double x,a,b,c; 

{ return{ a*b*c*pow(x.c-1)*exp(-b*pow{x,c)) ); } 

void sim_nhpp(lambda_nhpp. par1_nhpp, par2.nhpp, par3_nhpp, 

start.time. stop_time, path, nojx>ints) double lambda_nhpp{doubte,double,doubie,doubte); double 
^ par1_nhpp. par2_nhpp, par3_nhpp. start Jime, stop_time; double pathQ; long *no_potnts; 

double s(m_time,x,u,x_bar,lambda_bar; 
*no _points=0; 
sinnjlme = start_time; 

do{ 

if (par3_nhpp<=1) { lambda_bar = Iambda_nhpp(sim_time,par1_nhpp.par2_nhpp.par3_nhpp); 
else { 

x_bar=pow((par3_nhpp-1.0)/par2_nhpp/par3 nhpp,1,0/paf3 nhpp)* 
if (sim_time>=x_bar) { 

iambda.bar - Ianibda_nhpp(sfmjime,par1_nhpp.par2_nhpp,par3_nhpp); 
eise { 

iambda_bar = iambda.nhpp(x_bar,par1_nhpp,par2„nhpp,par3_nhpp); 
} 

X = sim_exp([ambda_bar); 
u = drand48{); 



if {u<=lambda_nhpp(sirn _time+x.par1_nhppppar2_nhpp.par3_nhpp)/IaiTibda_bar) { Cnojoints)**; 
pathCno_potnts]=sim time+x; 

} 

srm time+=x; 
} " 

while (sim_tinie<=stop_time); 

return; 

) 
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r00841 



r00851 



r00861 



Program 3: 

/* Definition of confidence ranges from tlie simulation data */ 
/* (tlie simulation data is sorted into ascending order) */ 
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#inciude <stdlib.h> 
#(ndude <math.h> 
#indude <stdio.h> 

int qsortJcmp(int*,int*); 
int qsortJcmp(x,y) 
int *x, *y; 
{ 

if (•x<*y) { return (-1 );} 
else rf (•x==*y) { return { 0 ); } 
else { return ( 1 ); } 

} 

int main(argc,argv) 
int argc; 
char *argv[]; 
{ 

mt pnt[11][1 00000]; 
int qs[1 00000]; 
char *dname; 
int fracj; 

long lautlower_bound,upper_bound; 
long l,no_pnt5eed_run: 
double ctime[11],x: 
FILE *datei; 



if (argcO) { 

printfC^nVnZuwenig Argumentel \n\n-); prlntfCAufruf: %s <Dateiname> <Konfiden2niveau (in 
%%)>\^\n^ argv[0]); return 1; 

} 

dname = argv[1]; 

frac = 100-atoi(argv[2]); 

lauf=0; 

date! = fopen(dname,"r); 
if (datei==NULL){ return 1;} 
e!se{ 

while {ifeof(datei)) { 
lauf++; 

for(i=1;i<=9;i++){ 

fscanfCdate^-yoSifyoSd". &ctime[i], &pntli][lauf]); 
fscanfCdatei/VoSIf %6d \n^ &ctinne[10], &pnt[10][lauf]); 

fcl05e(datei); 
} 

lower^bound = (long)fioor(!aurfrac/200.); upper_bound = {!ong)ceil(lauf*(200.-frac)/200 )• 
if (lower_bound<1 ) {lower_bound = 1 ;} ^ / /. 

r00871 
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printf("\n\n%2d%%-Sicherheitsbereich bei %d Simulationslaufen\n\n", 
lOO-frac,iauf); 
for (i=:1;i< ;!++){ 
for (l=1;l<=lauf;l++){ 
qsri} = pnt[i][l]; 
} 

qsort(&qsI13, lauf. sizeof(int), &qsortJcmp); 
printfC'Zeitpunkt: %8.2f Fehlerintervall: {%d,%d]\n". 
ctimefi], qs[lower_boundj, qs[upper_bound]); 



return 0 ; 
} 



12088] The above-described method and apparatus are illustrative of the 
principles of the present invention. [Referoncos: ] 

[P4 ] [Sidn e y I. Resnick: "Adv enturos][4f)-][ Stochastic Proc o ss o o", 

Birkhaus o r Booton, 1902, I SBN 3 76 1 3 3501 2, pp. 303 317] 

[P3 Bratoly ot a l ., 1087 ] 

Numerous modifi cations and adaptations will be readilv apparent to those 
skilled i n this art without departing from the spirit and scope of the present 

invention. 
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[Ab st ract ] 



[ Method and arrangement for predicting moasurement data using given 
measur e m e nt dota l ABSTRACT 

r00891 A method and an arrangement are [ proposod l provided for predicting 
5 measurement data using given measurement data, in which a stochastic process is 
matched to the given measurement data. Simulation runs are carried out from a 
given time-point until a final time-point. The forecast measurement data is 
determined for each simulation mn. Measurement data for the final time-point is 
predicted within a range of values, which is determined by the forecast 
1 0 measurement data. 
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Description 

Method and arrangement for predicting measurement data 
using given measurement data 

5 

The invention relates to a method and arrangement for 
predicting measurement data using given measurement 
data . 

I 

10 A technical system often requires facilities for 
forecasting based on known (measurement) data, 
particularly in the context of error susceptibility or 
cost estimates. 



15 Forecasts generated by experts are generally subject to 
errors. Experts cannot carry out exact analyses, at 
least of highly complex systems. 

A stochastic point process, in particular a Poisson 
20 process, is described in [1] . 

The object of the invention is to allow the automatic 
prediction (forecast) of measurement data using given 
measurement data. 

25 

This object is achieved in accordance with the features 
of the independent patent claims . Developments of the 
invention are described in the dependent claims. 

3 0 In order to achieve the object, a method is proposed 
for predicting measurement data using given measurement 
data, in which a stochastic process is matched to the 
given measurement data. Simulation runs are carried out 
from a given time-point until a final time-point. The 

3 5 forecast measurement data is determined for each 
simulation run. Measurement data for the final time- 
point 
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is predicted within a range of values, which is 
governed by the forecast measurement data. 

One development is to define a confidence range for the 
5 prediction of measurement data, where the a% lowest and 
b% highest forecast measurement data are eliminated. In 
particular, a% can equal b% . For example, a 95% 
confidence range can thus be defined by ignoring the 
2.5% lowest and 2.5% highest forecast measurement data. 

10 

One advantage is that the measurement data can be 
predicted (forecast) with an accuracy that is within a 
confidence range, from a given time-point. This makes 
it possible to identify e.g. the feasibility or 

15 impossibility of a task associated with the measurement 
data, at an early stage. Appropriate measures can 
therefore be initiated in order to counteract forecast 
impossibility. This is particularly important in the 
case of a complex system, e.g. a software development 

2 0 process, where the extent to which a schedule can be 
followed before the software is completed can be shown 
in a subsequent test phase. Even more important in this 
context is the ability to adopt countermeasures at an 
early stage if a delay has been clearly identified, 

25 e.g. in an integration test phase. This firstly affects 
the feasibility of the specified deadline (timescale) 
and secondly directly affects costs, since non- 
compliance with the agreed timescale often results in 
additional costs . 

30 

One refinement is for the stochastic process to be a 
non- homogeneous Poisson process. 

In particular, the measurement data may in one 
35 refinement comprise numbers of errors. This applies to 
software development , 
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for example, where the level of maturity is documented 
in accordance with the errors measured in a test phase. 
Completion is directly dependent on this level of 
maturity. In other words, the software cannot be 
5 delivered to customers until most of the errors have 
been removed from the software. This is particularly 
important with regard to resources (required to test 
and correct errors) and costs (due to delayed 
delivery) . 

10 

In order to achieve the object, a method is also 
proposed for predicting measurement data using given 
measurement data, in which a stochastic process is 
matched to the given measurement data. A range is 
15 ascertained, by sorting the probability values 
generated by the stochastic process according to size, 
around an expected value. Measurement data is predicted 
on the basis of this range, and in particular the 
probability values within the range. 

20 

One development is for the probability values generated 
by the stochastic process to be sorted symmetrically by 
size around the expected value. In particular, this 
means that the highest probability value represents the 

2 5 middle of the range, i.e. the expected value, whereas 

the next highest probability value is arranged to the 
right or left of the expected value. The next highest 
probability value is then arranged symmetrically on the 
other side of the expected value, in turn. 

30 

This analytical (design) procedure provides a range, 
where the breadth of the range in turn indicates which 
probability values are significant in the prediction of 
the measurement data . 

3 5 

In one particular refinement, the breadth of the range 
is determined by ignoring the 
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probability values that lie below a given threshold. 

This produces a range (confidence range) , which has a 
specific breadth as a result of the threshold. This 
breadth corresponds to the certainty with which the 
measurement data is predicted. 

If one assumes that the stochastic process is a non- 
homogeneous Poisson process, then the non- homogeneous 
Poisson process defines a step size, particularly on a 
time axis t, which indicates when the next error will 
occur. One characteristic of the non- homogeneous 
Poisson process is that it has no memory, so that a 
"no-memory" search is carried out from each error that 
occurs at a specific time-point, for a time-point that 
indicates the next error. 

In order to achieve the object, an arrangement is also 
proposed for predicting measurement data using given 
measurement data, whereby a processor unit is provided 
and configured in such a way that: 

a) a stochastic process can be matched to the given 
measurement data; 

b) simulation runs can be carried out from a given time- 
point until a final time-point; 

c) the forecast measurement data can be determined for 
each simulation run; 

d) the prediction of measurement data for the final 
time-point can be predicted within a range of values, 
which is determined by the forecast measurement data. 

In order to achieve the object, an arrangement is 
further proposed for predicting measurement data using 
given measurement data, whereby a processor unit is 
provided and configured in such a way that: 
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a) a stochastic process can be matched to the given 
measurement data; 

b) a range can be ascertained by sorting probability 
values generated by the stochastic process according 

5 to size around an expected value; 

c) the measurement data is predicted within the limits 
of the range . 



The arrangements are particularly suitable for carrying 
10 out the inventive method or the developments described 
above . 



Exemplary embodiments of the invention are shown and 
explained below with reference to the drawings, in 
15 which: 



Fig, 1 is a graph showing an accumulated number of 
errors over a test period; 

2 0 Fig. 2 is a graph showing the superimposed confidence 
ranges for different process models; 

Fig, 3 is a flowchart showing the steps in a method 
for predicting measurement data using given 
2 5 measurement data ; 



30 



Fig. 4 is a further flowchart showing the steps in a 
method for predicting measurement data using 
given measurement data; 

Fig. 5 shows a processor unit. 



35 



In order to be able to forecast a number of errors to 
be expected in a technical process, e.g. in a software 
development process, non-homogeneous Poisson 
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processes (NHPP) are 
measurement data, e.g. 
time, as follows: 



- 6 - 

calibrated, i.e. matched to 
the occurrence of errors over 



5 The following equation describes a counting process 
associated with the stochastic point process (non- 
homogeneous Poisson process) : 

{N(t)}tgR+ (1) 

and a time-point to defines the end of a test period, 
10 i.e. a time -point at which the given data ends. The 
stochastic processes 

W}t6R+ and (2) 
{0{t)}^^jj+ (3) 

are searched with 
p(u(t) < N(t) - N(to) < 0(t) I N(to) = no) > a (4) , 

15 for all time-points where t>to and given values cc g 
(0.1) (confidence level) and no e N. In particular, the 
following text examines the increases in the stochastic 
countings process in relation to the time-point to- 

2 0 In the present case, where equation (1) represents a 
non- homogeneous Poisson process, the following equation 
(cf. [1]) 



25 



Mtl) - 4t0) = <) = exp(- [i(ti) - i(to)]) • 



(5) 



applies for 

0 < to < ti < CO, ^ € No (6) 

and an intensity (mean measure, mean value function) of 
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i: R"^ R"^,t i-^ i{t>= EN(t) 



(7) 



Since the nature of the Poisson process dictates that 
the increases (error increases in this case) are 
independent of previous increases, equation (5) for the 
time-points t>to to define a (minimum) range 



[gu' go] = [gu(t), go(t)] c= Nq 

can be simplified to 

go 

2 P(N(t) - N(to) = e) > a 

^=gu 



(8) 



(9) . 



Due to the unimodal nature of the Poisson count 
density, a range [gu, go] can be determined as follows: 



Step 1: Sort the elementairy probabilities 

: = p(N(t) - N(to) = i)r i G No 

into descending order and label the values 
sorted thus using 

P(0)'P(1)'--- {PO' PI' • • •} = {p(0)' P(l)' • • •) 

P(0) ^ P(i) 2: . . . ) ; 



step 2: Determine ^min ■ " inim^ € Nq 



Z p(i) ^ « 

1 = 0 



Step 3 : Determine an index set 



: = jio, . . . , Ur^i^ } No where 
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min 



Step 4: Substitute := min{i} and 



The range from equation (8) is also referred to as the 
forecast range . 

Stochastic simulation (second approach) 

It is possible to determine the confidence range 
described using simulation, with the following steps: 

Step 1: Start independent simulation runs based on the 
selected process model at time-point to of the last 
error message m £ N; 

Step 2: End a simulation run as soon as the required 
final time-point te is reached; 

Step 3 : Repeat Step 2 until all simulation runs are 
finished; 

Step 4: Sort the numbers Ni(te) of the errors generated 
in the i-th simulation run in the time period (to, te) / 
i = l, m, in descending order, and label the values 

sorted thus ^(-^^(t^), . , . , N(^)(tg) ; 

Step 5: Substitute 




and 



GR 98 P 5929 

- 9 - 

i.e. eliminate the (100? (1-a) /2) % lowest and 
highest values. 

This produces the confidence range directly. 

Each individual simulation run is based on a simulation 
algorithm, which is known from (cf. [2]) : 

Simulated generation of intermediate arrival times for 
a non- homogeneous Poisson process: 

Step 1: Substitute X sup{A.(t)}, where: 

t>tc. 



(10) 



step 2 : Generate a (pseudo) random variable X that is 
exponentially distributed with the parameter X, i.e. 
X := - log(u) / X , where U is equally distributed over 

(0,1) . 

Step 3 : Generate a random variable U that is equally 
distributed over (0,1) . 

Step 4: If U < X{t^ + x)/X, then substitute t* : = ts + X; 
otherwise substitute ts := tg + X and go to Step 1. 

The example graph in Fig, 1 shows an accumulated number 
of errors during a given test period. From time-point 
to. it shows a prediction range for all time-points to + 

X. 



The intensity i is normally derived from equation (10) 
for X. For example the result is as follows: 
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a) X(t) ^ a • b • c • exp(- bt^) • t*^"^ 

{X{t) is Strictly monotonously descending for 
c < 1, and unimodal for c > 1 with a definitive 
maximum at a point 
/c - 1 , 



' max 1/ 

be 



b) Otherwise, X is derived in accordance with the 
above comments as follows: 

- ^ M^sl (c ^ 1) V (ts ^ tmax) 



The graph in Fig. 2 shows the superimposed confidence 
ranges. In particular, this illustrates that possible 
forecasts become more scattered the further they extend 
into the future. In particular, confidence ranges 
calculated using different process models can be 
demonstrated in the same way as shown in Fig. 2. 

Fig. 3 shows a flowchart for the steps of a method for 
predicting measurement data using given measurement 
data. In Step 301, a stochastic process, in particular 
a non-homogenous Poisson process (to represent a 
stochastic count process) , is matched to given 
measurement data. In Step 3 02, simulation runs are run 
from time-point to to a final time-point te that is to 
be forecast. In Step 303, for each simulation run, 
forecast measurement data is determined and a 
prediction of measurement data is restricted to a range 
which is covered by the measurement data determined by 
the simulation runs (see Step 304) . In Step 305, a 
confidence range is determined in which a given 
proportion of the lowest and highest forecast 
measurement data is ignored in each case (this 
corresponds to the aforementioned range) . The method 
terminates in Step 3 06. 



GR 98 P 5929 

- 11 - 

Fig, 4 shows a further flowchart for the steps of a 
method for predicting measurement data using given 
measurement data. In Step 401, a stochastic process, in 
particular a non- homogenous Poisson process, is matched 
to the given measurement data. Probability values are 
determined using the stochastic process , and these are 
sorted according to size around an expected value (see 
Step 402) . This sort operation results in the 
definition of a range, namely a confidence range in 
this case. The breadth of the confidence range is 
determined by comparing the accumulated probabilities 
with a given threshold. As described above, the 
confidence range gives a distribution or uncertainty, 
respectively, of a time-point to in the future, which 
allows the measurement data to be estimated in the 
future (see Step 403) . The method terminates in Step 
404 . 

Fig. 5 shows a processor unit PRZE . The processor unit 
PRZE comprises a processor CPU, a memory unit MEM, and 
an input/output interface lOS, which is used in 
different ways via an interface IPC: a graphics 
interface allows output to be viewed on a monitor MON 
and/or output to a printer PRT. Inputs are entered via 
a mouse MAS or a keyboard TAST. The processor unit PRZE 
also includes a data bus BUS, which provides the 
connection between a memory unit MEM, the processor CPU 
and the input/output interface lOS. It is also possible 
to connect additional components to the data bus BUS, 
e.g. additional memory, data storage (hard disk) or 
scanner . 

The C programming language is used in the following 
examples, which show an algorithm to define confidence 
ranges for forecasts and an algorithm for simulated 
definition of confidence ranges for forecasts. 
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/* Definition of confidence ranges for forecasts */ 

/* based on the generalized Goel-Okomoto model */ 

#incliide <stdlib.h> 
#(nclude <rnath.h> 
#include <stdio.h> 

#define true 1 
#define false -1 

double mv_genGO(doubie,doubie.doub!e,double); double poisson(doublejong): void ki^nhppO 

int main(argc,argv) 

int argc; 

char *argv[}; 

{ 

double a.b,c,bt,st,kn; 
long low.upp.Jauf; 

if (argc<7) { 

pr[ntf("\n\nZuwer!ig Argumente! \n\n"); 

printffAufruf; %s <Par1> <Paf2> <Par3> <Startzeit> <End2eit>", 
"<KNiveau>\n\n", argv[0]); return 1; 
} 

a = atof(argv[1]); 
b = atof(argv[2]); 
c = atof{argv[3]); 
bt= atof(argv[4]); 
st= atof(argvI5]); 
kn= atof{argv[6]); 

for (lauf=1;lauf< ;iauf++) { 

kLnhpp(mv_genGO.a,b,c,btbt+laur(st-bt)/10.,kn,&low,aiupp); 
printf("ZeJtpunkt: %8.2f Fehlerintervall: [%d,%d]\n". 
bt+laur(st-bt)/10.. low, upp): 

} 

return 0; 
} 

double mv_genGO(x,a,b,c) 
double x.a,b,c; 

{ retum( a*{1.0-exp{-b*pow(x,c))) ); } 

double poisson{iambda,wert) 
double lambda; 
long wert; 
{ 

long i; 

double itval.hv; 

if (lambda<600) { 
itval = exp(-lambda); 

for (i=wert;i>=1;i-) { itval *= lambda/(double)i; ) 
} 
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ei$e{ 

hv = exp(-lambda/(double)wert); 
itvat = 1.0; 

for (i5=wert;i>=1;i-) { rtvai *= lambda/{double)i*hv; } 
} 

return ( rtval ); 
} 

void kLnhpp(mv_nhpp, par1_nhpp, par2„nhpp, par3_nhpp. 

start_trme. stop^time, k_niveau, lower, upper) double mv_nhpp(doubie,double,doubie,doubie): 
double par1_nhpp, par2 nhpp. par3 nhpp» start time, stop time, k_niveau: long *lower, *upper 

{ 

long iauf; 

int lborder,modJow,mod_upp; 
double sum,tmp_mv, valj, vat_u; 

tmp_mv - mv_nhpp{stop_time,par1_nhpp,par2_nhpp,par3_nhpp) - 
mv_nhpp(start_time,par1_nhpp,par2_nhpp,par3_nhpp); lauf = (long)tmp_mv; 
*low€r = iauf-1 ; 

*upper = iauf+1; modJow= false; mod„upp= false; sum = poisson(tmp„mv,lauf); valJ = 

poisson(tmp_mv,*iower); val_u = polsson(tmp_mv *upper); 
while (sum<k_niveau) { 
if {valj<val_u) { 
sum += val_u; 
{*upper)++; 
Iborder = false; 
mod_upp = true; 

val_u = poisson(tmp_mv,*upper); 
} 

else { 

sum valJ; 
(*iower)— ; 
Iborder - true; 
mod Jow = true; 

val I = poisson(tmp mv *lower); 
f 

} 

if (Iborder == true) { (•|ower)++: } 
else { Cupper)-; } 

if (modjow == false) { (*lower)++; } 
if (mod_upp false) { (*upper)-; } 

return; 
} 
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/* Simulated definition of confidence ranges for 
forecasts */ 

/* based on the generalized Goel-Okomoto model */ 

#include <stcilib.h> 
#!nciude <math.h> 
#include <time.h> 
#inciude <stdio.h> 
#inc{ude <va)ues.h> 

#define taie 1 
#define false -1 

double drand48(votd); 
void srand48(long); 

double sim_exp(double); double lambda_genGO{double,double,double,double); void sim_nhpp(); 

int main(argc,argv) 

int argc; 

char 'argvQ; 

{ 

time_tt; double a,b,c,btst,pnt[1000000].check time[12]; long lauf»no ^pnt.seed run; intclauf; 
FILE *datei; 
if (argc<6) { 

printf{'^n\n2uwenig Argumente! \n\n"); 

printffAufnjf. %s <Par1> <Par2> <Par3> <Startzeit> <Endzeit>\n\n", 
argvIO]); return 1 ; 
} 

datei = fopenC'sim.seed\"r"); 
if {datei==NULL) { 
seed_run = 1 ; 

} 

else { 

fscanf(datei,"%6d",&seed_run); fclose(datei); seed_run++; 
} 

datei = fopenrsim.seed*,''w+''); 
fprintf(datei, "%6d\n", seed_run ); 
fclose(datei}; 

time (&t) ; r initiatisierung des */ 

t += seed_run*100 ; /* Zufallszalilengenerators V 
srand48 {(unsigned long) t) ; /* mit Hilfe der Systemzeit */ 
a = atof(argv[1]); 
b atof(argv[2]); 
c = atof(argv[3]); 
bt= atof(argv[4]); 
st= atof(argv[51); 

stm_nhpp(lambda_genGO,a,b,c,bt,st&pnl,&no_j3nt); 
for (lauf=1 ;iauf<=no _pnt;Jauf++) { 
pnntf("%15.7f %10d \n^ pnt[iaufl, lauf); 
} 
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datei - fopen("ki.tmp",*a''); 
foT{lauf=1:iauf<;lauf++){ 

check_time[laaf] = bt+laLtf*(st-bt)/10.: 

check_t!me[11] = pntlno_pnt]+1; /* groaer ais die groate 
simulierte Zeit */ 
ciauf - 1; 

for (lauf=1;lauf<=no _pnt;lauf++) { 

while (pnt[lauf)>=check_time[claufl) { fprintf(datei, "%8.2f %6d checkjimelclauf], tauf-1); clauf+-»-: 
} 

} 

if (pnt[no_pnt] <check timeJIO]) { 

for (lauf=clauf:lauf< ;lauf++) { fprintf(datei, "%8.2f %6d checkjimeliauf], no_pnt); 
} 

} 

fprintf{datei, "\n''); 
fclose(datet); 

return 0; 
} 

double sim_exp(lambda) 
double lambda; 

{ r€turn( -tog{drand48{))/iambda ); } 

double lambda jgenGO(x,a.b,c) 
double x,a,b,c; 

{ retum( a*b*c*pow(x,c-1)*exp(-b*pow(x,c)) ); } 

void sim_nhpp(lambda_nhpp, par1_nhpp, par2_nhpp, par3_nhpp, 

staft_time. stop_time, path, no_points) double lambda_nhpp(double, double. double, double); double 
par1_nhpp, paf2_nhpp, par3_nhpp, startjime, stopjime; double pathQ; long •no_points; 

{ 

double sim_time,x,u,x_bar,lambda_bar; 

*no_points=0; 

sim_time = start^time: 

do{ 

If (par3_nhpp<=1) { lambda_bar = Iambda_nhpp(simjime,par1_nhpp,par2_nhpp.par3_nhpp); 
} 

else { 

x_bar =: pow({par3_nhpp-1.0)/par2_nhpp/par3_nhpp,1.0ypar3_nhpp); 
if (sim_tim6>-x_bar) { 

lambda bar « lambda„nhpp(sim tlme,par1_nhpp,par2 nhpp,par3_nhpp); 
} 

else { 

lambda bar = lambda nhpp(x bar,par1 nhpp,par2 nhpp,par3_nhpp); 
} " 

} 

x = sim_exp(lambda_bar); 
u = drand48(); 
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if {u<=lambda_nhpp{5im_time+x,par1_nhpp,par2„nhpp,par3_nhpp)/lamfada_.bar) { (*no_points)' 
path[*no_points]=sim time+x; 

} 

simjme-^^x; 
} 

while (sim_time<=5stop_time); 

return; 

} 
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/* Definition of confidence ranges from the simulation 
data */ 

/* (the simulation data is sorted into ascending order) 

#tnclude <stdlib.h> 
#inciude <math.h> 
#include <stdio.h> 

int qsortJcmp(inf,int*); 
int qsortJcmp(x,y) 
int *x, *y; 

{ 

if (*x<*y) { return ( -1 ); } 
else if (*x==*y) { return { 0 ); } 
else { return ( 1 ); } 

} 

int main(argc,argv) 
int argc; 
char *argv[]; 

{ 

int pnt[11][1 00000]; 
intqs(100000]; 
char *dname; 
int fraci; 

long tauf,lower_bound,upper_bound; 
long l,no_pnt,seed_run; 
double ctime[11],x; 
FILE -datei; 



If (argcO) { 

pnntfC'\n\nZuwenig Argumentel \r\\ny, printf("Aufruf: %s <Dateiname> <Konfiden2niveau {in 
%%)>\n\n-, argv[01); return 1; 

} 

dname = argv[1l; 

frac = 100-atoi(argv(2]); 

lauf = 0; 

datei = fopen(dname, V); 
if (datei==NULL) { return 1; } 
else{ 

while (!feof(datei)) { 
iauf+-*-; 

for(i=1;i<=9;t++){ 

fscanf(datei;%8if%6d". &ctinr?e[i]. &pnt[i][lauf]): 
} 

fscanf(dater;%8lf %6d \n-, &ctime[10], &pnt[10][lauf]): 
} 

fclose(datei); 
} 

lower_bound = (long)floor{laurfrac/200.); upper_bound = (Iong)ceil(laun200,-frac)/200.); 
if (lower_bound<1) {lower_bound = 1;} 
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printf("\n\n%2d%%-S}chertieitsbereich bei %d Simulationslaiufen\n\n". 

100-frac,lauf); 

for (i=1;i< 

for (l=1;l<=lauf;l++){ 
qs[n = pnt[il[l]; 

} 

qsort(&qsI1], lauf, sizeof(lnt), &qsortJcmp); 
printfC'Zeitpunkt: %8.2f Fehlerintervall: {%d,%dl\n-. 
ctime[l], qs[Iower_bound], qs[upper_bound]); 

} 



return 0 ; 
} 
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1. A method for predicting measurement data until a 
final time-point using given measurement data, in 
which: 

a) a stochastic process is matched to the given 
measurement data; 

b) simulation runs of the stochastic process are 
carried out from a given time-point until the 
final time-point; 

c) the forecast measurement data is determined for 
each simulation run; 

d) measurement data is predicted by stating a range 
of values, which is determined by the forecast 
measurement data . 

2. The method as claimed in claim 1, in which: 

a confidence range is determined for the 
prediction of measurement data, where the a% 
lowest and b% highest forecast measurement data 
are eliminated . 

3. The method as claimed in claim 2, in which: 
a% and b% are equal . 

4 . The method as claimed in one of the preceding 
claims, in which: 

the stochastic process is a non- homogeneous 
Poisson process. 

5. The method as claimed in one of the preceding 
claims, in which: 

the measurement data represents numbers of errors. 

6. The method for predicting measurement data using 
given measurement data, in which: 

a) a stochastic process is matched to the given 
measurement data; 
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b) a range is ascertained, by sorting the probability- 
values generated by the 
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stochastic process according to size, around an 
expected value; 
c) measurement data is predicted within the limits of 
the range . 

The method as claimed in claim 6, in which: 
the probability values generated by the stochastic 
process are sorted symmetrically by size around 
the expected value . 

An arrangement for predicting measurement data 
until a final time-point using given measurement 
data, whereby a processor unit is provided and 
configured in such a way that : 

a) a stochastic process can be matched to the given 
measurement data ; 

b) simulation runs of the stochastic process can be 
carried out from a given time-point until the 
final time -point ; 

c) the forecast measurement data can be determined 
for each simulation run; 

d) measurement data is predicted by stating a range 
of values, which is determined by the forecast 
measurement data . 

An arrangement for predicting measurement data 
using given measurement data, whereby a processor 
unit is provided and configured in such a way 
that: 

a) a stochastic process can be matched to the given 
measurement data; 

b) a range can be ascertained by sorting probability 
values generated by the stochastic process 
according to size around an expected value; 
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c) the measurement data is predicted within the 
limits of the range. 
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Abstract 

Method and arrangement for predicting measurement data 
using given measurement data 

A method and an arrangement are proposed for predicting 
measurement data using given measurement data, in which 
a stochastic process is matched to the given 
measurement data. Simulation runs are carried out from 
a given time-point until a final time-point. The 
forecast measurement data is determined for each 
simulation run. Measurement data for the final time- 
point is predicted within a range of values, which is 
determined by the forecast measurement data. 
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VERTRETUNGSVOLLMACHT: Als benannter Erfinder 
beauftrage ich hiermit den nachstehend benannten 
Patentanwalt (oder die nachstehend benannten 
Patentanwalte) und/oder Patent-Agenten mit der 
Verfolgung der vorliegenden Patentanmeldung sowie 
mit der Abwicklung aller damit verbundenen Geschafte 
vor dem Patent- und Warenzeichenamt: (Name and 
Registrationsnummer anfuhren) 



German Language Declaration 



POWER OF ATTORNEY: As a nanned inventor, 1 
hereby appoint the following attorney(s) and/or 
agent(s) to prosecute this application and transact all 
business in the Patent and Trademark Office 
connected therewith, (list name and registration 
number) 



And I hereby appoint 



Customer No. 26574 



Teiefongesprache bitte richten an: 
(Name und Telefonnummer) 



Direct Telephone Calls to: (name and telephone 
number) 



Ext. 



Postanschrift: 



Send Correspondence to: 

Schi_f f.JHaLdin& Waite 
6600 Sears Tower 60606-6473Jau[aa£ia Illinois 
Telephone: +131^25§_5ZSa^cl Facsimile +1 312 258 5921^ 



Custo 



No. 26574 



Voller Name des einzigen oder ursprunglichen Erfinders. 

s THOMAS ACKERMANN. 


Full name of sole or first inventor: 

XHQMAS ACKFEMANM 


UnterscbrlfrnSf Erfinders Datum 


Inventor's signature Date 


Wohnsitz ' I/' 

HERZ0J3O4AL1SACM. DEUTSCHLAND 


Residence 

:|-|FP7nnFMAI IRAP.H GFRMANY -^^^^ 


Slaatsangehorigkeit 

DE 


Citizenship 

DE ^ 


Postanschrift 

SCHULSTRASSE 28 


Post Office Addess 

SCHULSTRASSE 28 


91074 HERZOGENAURACH 


91074 HERZOGENAURACH 


Voller Name des zweiten Miterfinders (falls zutreffend): 

Dr. MICHAEL GREIMER 


Full name of second joint inventor, if any: 

Dr. MICHAEL GREINER 


Unterschrift des Erfinders Datum 


Second Inventor's signature Date 


Wohnsitz ^ 

DACHALLDEUTSCHLAND 


Residence 

DAChlAU, GERMANY £Pt^ 


StaatsangehOrigkeit 
DE ^ 


Citizenship 

DE ^ 


Postanschrift 

PRAELAT-WOLKER-STR. 2 


Post Office Address 

PRAELAT-WOLKER-STR. 2 


85221 DACHAU 


85221 DACHAU 



{Bine entsprechende Informationen und Unterschhften im 
Falle von dritten und weiteren Miterfindern angeben). 



subsequent joint inventors). 
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Voller Name d&s dritten Miterfinders! 

Dr. PETER LIGGESMEYEE. 


Full name of third joint inventor: 

Dr. PETER LIGGESMEYER 


Untetsi^rift ctes Erfinders Datum 


Inventor's signature 


Date 


Wohnsitz fC^'^^ O*^ ( 

Potsdam, Deutschland 


Residence 

Potsdam. Germany £> 


Staats angehoriQkeit 

DE ^ 


Citizenship 

DE 


Postanschrift 

Eichenring 27 


Post Office Address 

Eichenring 27 




14469 Potsdam 


Voller Nam© des vierten Miterfindersi 

OLIVER iVlAECKELc 


Full name of fourth joint inventor: 

OLIVER MAECKEL 


1 Inttfircrhrift riA<i Frfinrlpr<5 .^^''^^^^O Datum 


Inventor's signature 


Date 


MUENdHElSl, DEUTSCHLAND & <PX 


Residence 

MUFNCHEN^GERMANY 


OlacLXoanycflUi lyKcIL 

DE*^ 


Citizenship 

DE. 


Postanschrift 

SALMDORFER STR. 5 


Post Office Address 

SALMDORFER STR. 5 


81929 MUENCHEN 


81929 MUENCHEN 


VriHpr Mamp rip« funftpn Mi+prfinders' 

VUlJwl INdlllv? UCO lUIIILwII IVI J Lwl 1 M l\4wt 


Full name of fifth joint inventor: 


unierscnnn u€b cninutjio uaium 


Inventor's signature 


Date 


Wohnsitz 


Residence 

3 




Citizenship 


Postanschrift 


Post Office Address 






Voller Nsme des sechsten Miterfindersi 


Full name of sixth joint inventor: 


1 Intdrcr-hriff Ha<; FrfinriAr<! Datum 


Inventor's signature 


Date 


VV \j\ 11 lOlu^ 

I 


Residence 




Citizenship 


Postanschrift 


Post Office Address 







{Bitte entsprechende Informationen and Unterschriften im (Supply similar information and signature for third and 
Falle von dritten und weiteren Miterfindern angeben). subsequent Joint inventors). 
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